Asymptotic analysis of the model for distribution of high-tax payers 
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The z-transform technique is used to investigate the model for distribution of high-tax payers, 
which is proposed by two of the authors (K. Y and S. M) and others|6,, 7,, 8]. Our analysis shows 
an asymptotic power-law of this model with the exponent -5/2 when a total "mass" has a certain 
critical value. Below the critical value, the system exhibits an ordinary critical behavior, and scaling 
relations hold. Above the threshold, numerical simulations show that a power-law distribution coex- 
. ists with a huge "monopolized" member. It is argued that these behaviors are observed universally 

' in conserved aggregation processes, by analizing an extended model. 
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O 1. Introduction 

Some years ago, two of the authors (K. Y and S. M) et al [1, 0, S] proposed a simple model for distribution of 
P\J ] high-tax payers. They considered a 2-body interaction system where the winner takes all competing money when the 
^ ' competition occurs. Their model is regarded as a simple aggregation system, which is widely discovered in various 
, systems in nature and, naturally, has been attracted considerable interest by many researchers 0,0] • Their simulation 

■ gives the universal power-law distribution, and shows a good fit with the date in Japan and the United States' CEO's. 
, However they could not succeed in theoretical derivation of obtained results. 

7-H ■ We analyze their model mathematically, and derive an asymptotic expression of the probability distribution func- 
i tions. In the model, the total number of units (mass) is conserved, and plays a role of a control parameter. The 

■ distribution function obeys the power-law at a certain critical value. Below the critical value, the system exhibits an 
ordinary critical behavior, and scaling relations hold. Above the threshold, numerical simulations show that a power- 
law distribution coexists with a huge "monopolized" member. These bahaviors are agreed with the phase transision 

G ' of the non-conserved process discussed by Krapivsky and Render[l| and Majumdar et al^. 

O 2. Model and Analysis 

The model introduced by Yamamoto et al, is defined by the following competition processes and conservation 
conditions: 



CI. The system consists of N homogeneous members. The number N is fixed. 

C2. The system holds S units as resource money. Each member holds one unit as a minimum amount. The total 
amount S is conserved. 

PI. Two members, who are picked at random from the whole membership transact an economic activity competi- 
tively. The winner collects the whole of the competing money and the loser loses all one's money. In order to 
keep the number of active members constant, one unit is added to the loser. 

P2. To preserve the total resource money, one unit is reduced from a member who is picked at random and has 
resource money more than or equal to 2 units. 
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We analyze this model by using the z-transform technique and calculate the probability distribution function P{X) 
that an arbitrarily chosen member holds X units of resource money. At steady-states, the equations for P{X) read 



(1-F(1))(1-P(1)) + 



P{2) 



={pii)r (x = 1), 



l-p(l) 



(1) 

(2) 



Note that J2x ^(^) — 1 ^^r normalization. For the analysis of the equation of this kind, it is useful to consider the 
z-transform (j){z) of P{X): 



(Z) = £ P{X)2 



-X 



(3) 



Here it is defined that P(0) = 0. From the normalization condition, 0(z) must satisfy 0(1) = 1. Taking the z-transform 
of Eq.(l2|) and Eq.([T]), we get the basic equations in z-space as, 



{1 - P{l)}(jP{z) + {z - 3 + 2P(l)}0(z) - P(l) + - = 0. 



(4) 



We can now analyze P(l) by differentiation Eq.((3]) with respect to z. Differentiating once and substituting z = 1, 
we simply get the identity. Next differentiation leads to the equation about 0' = (/"'(l), 



{l-P(l)}0'2 + 0' + l = O. 



(5) 



Note that — 0'(1) = ^ XP{X) = S/N is the average of resource units, and is conserved in processes. We employ —(p' 
as the control parameter of the system. Once —0' is given, P(l) can be determined via Eq.(l5]). Fig. [T] shows P(l) 
versus —0'. It becomes evident that 1 > P(l) > |, and generally P(l) corresponds to two values (j>A and (j)B of —cj)', 
in the region A: 1 < (j>A < 2 and B: 2 < 0^. Note that from the definition of the system, each member should have 




FIG. 1: Value of P(l) versus the control parameter —0'. P(l) ~ 1 means that all the member have one unit. 



at least one unit thus —0' = 1 is the minimum value. 

First, we treat the case where —0' = 2 and Eq.([5|) has a repeated root, i.e., P(l) 
to ([4|, we have the solution of the quadratic equation 



3-2z±2z(l- 

z 



3/4. Substituting P(l) = 3/4 



(6) 



The singularity of 4>{z) near z ^ 1 is. 



0(z)U^i~l±2(l--)t 
z 



(7) 
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Then the asymptotic behavior of P{X) at large X is given by 

PW-^, (8) 

This power-law exponent coincides with numerical results of Q . 

Next we discuss the behavior of the system with the control parameter —0' 7^ 2. The original equations ([T|), ([2]) 
can be solved sequentially if once -P(l) is given. In Fig. [2l numerical solutions of Eq.((T]) and ([2]) are plotted for several 
values of P{1). With the values of — in the region A of Fig. [1] our simulation results agree with these sequential 



P(X) 




10° 10' 10^ 10^ 10'' 10' 10'^, 



FIG. 2: Solutions of P{X) with different values of P(l). (a) P(l) = 0.76, (b) P(l) = 0.751, (c) P(l) = 0.7501, (d) P(l) = 
0.75001, (e) P(l) = 0.75. 



solutions. We conclude that in this region the system has stable steady states, and the solutions of basic equations 
([D) and ([2]) represent these stable states. 

The behavior of these solutions, which is shown in Fig. [21 is very similar to the ones of ordinary critical phenomena. 
For arguing this behavior, let P{1) = | + e, then analyze the reverse z-transform. 



P{X) = ^ ^ 



(9) 



in the limit e — )■ 0. The most singular part of ^ is given by. 



p(X) - — / e-^«^-i)-^^+- >(z - \Wz-l + Aedz 
2Tn J 

^ J- i e^'^eaVea + 4e • eda (^^^ 
2m J 



- X-^^iieX), 
where the scaling function '^'{Y) is expressed as, 



*(^) = i> e°'^aVaT4da. (11) 

27^^ J 

We find that ordinary scaling relations hold. 

The behavior of P{X) is completely different when the control parameter —0' is larger than 2. Fig. [3] shows the 
simulation result for N — 100000 and — (/>' = 202. The part of units, about 2 x N units, obey the power-law with the 
exponent —5/2, and the rest of them coagulates eventually to one "monopolized" member, represented by the arrow 
in Fig. O Thus this member has almost ttir = {—(p' — 2) x N units. We have also done simulations with different 
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P(X) (timos 202N) 




10° 10^ 10^ lo' lo' 



FIG. 3: Simulation result for N = 100000 and -0' = 202. 



values of — in the range from 3 to 1002. In all these simulations P{X) behaves like Figs.[3l 

The behavior described above is substantially the same as that of the model reported by Majumdar et a/[^, though 
thier model is different from the present model for it does not preserve the total mass strictly. 



3. Extended Model 

Now we argue that systems which exhibit those behaviors are not limited to this particular process. To the end 
,we examine the model extended as follows, (i) In the competition process PI, 7i units are added to the loser, (ii) In 
the procedure P2, m member with X > 2 are picked at random and one unit is reduced from each of these member. 
The case n = m = 1 corresponds to the original model. 

Similarly to the original model, the equations for P{X) at steady state read 

""^^^^ =2F(1) (X = l), (12) 



l-P(l) 



P(^)P(J) + (l-P(n))(l-P(n)) + ^^— -^ = (P(n))2 + ^-A^ {X ^ n), (13) 

i+j=n ^ ' ^ ' 

E nmj) + ^^^^ - 2P(X) + iX^n,X>2). (14) 

i+3=X ^ ' ' 

(15) 

In z-space, these equations are rewritten as 

A(j)'^{z) + {m{z - 1) - 2A}(t){z) m{l - A) (- - l \ + — = (16) 

\ z I z" 



and the solution of ITBl) is 



1 



D{z) = {m{z - 1) - 2A}2 - 4A |m(l - A) Q - 1^ + ^| . (18) 

First we treat the case n ^ m. In this case the total amount S of units varies, thus the situation is quite different 
from the original model. However, the analysis described above is still applicable. 
When n m, the discriminant D{z) is 



D{z)^{z-l)f{z), /(1)^0. 



(19) 
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The singularity of ^(z) near z '-^ 1 is 0(^)12^1 '--^ 1 ± {(n — m){z — 1)}^/^. It gives a familiar power-law distribution 
(for example, see [3,|5|), 

P(X)^^, a=| (20) 

When n < m, the total amount S eventually vanishes. However if we take S large enough, our simulation shows 
the temporal power-law distribution with the exponent —3/2. Also in the case n > m, the range of the power-law 
distribution with the exponent —3/2 continually extends like y/i as the time t goes on. 

In the case n = m, — 0'(1) = —0' can be used as the control parameter, as the original model, since S is conserved. 
When -0' = 2, 

P{1) = 1- A = ^ ^ i2n-y/2n^-2n + i], (21) 
Z[n — 1) L J 

and 

D{z)^{z-lfg{z), g{l)^0. (22) 
Thus, the asymptotic behavior of P{X) at large X is given by 

P(X)~^, a^l- (23) 

It is just the same exponent as the original model. Also our simulations show that the behavior when —(p' ^ 2 is 
essentially the same as the original. 



4. Discussion 

These models show drastic change of behavior as the control parameter —0' varies. When —0' > 2, once the 
"monopolized" member appears our simulation shows this member does not exchange with other members after that. 
It is practically stable, so we believe that this transition is not only the normal static phase-transition, but also the 
dynamic "ergodic-nonergodic" transition. 

Our analysis described above shows that the exponent of "power-law" part does not change in two models. Further- 
more, it can be shown analytically that the exponent of the power-law depends only on conservation conditions when 
the kernel is constant. These studies suggest strongly that the behavior described above is observed universally in 
conserved aggregation processes. We believe this analysis gives clues for classification of these conserved aggregation 
processes. 

We thank Prof. S. Redner for pointing out references |lj and Q to us. 
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